Wandering in Color-Space* 
— why the life of pentaquark is so long ? — 
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The problem of the long life time of the pentaquark 0^ is investigated on the basis of the color 
molecular dynamics simulation. We find that it takes a long time (typically of 50 — 100 fm/c) for 
the initial pentaquark-state to rearrange its color and spatial positions to decay into the nucleon + 
kaon final state. Structure of the potential surface in the color and position spaces also supports 
this picture. Pentaquark wanders on the potential surface to find a narrow channel to decay. 
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^ : I. INTRODUCTION 

, An exotic baryon, pentaquark 0+, has been first reported by LEPS collaboration at SPring-S Later, there are 
other experiments reporting the confirmation of its existence, which is summarized in refs. [^|^. Since the 0"*" decays 
into a neutron (n), and a kaon (if^), its quark content is considered to be uudds. The mass is found to be about 
1540 MeV, while there is yet no direct measurement of its spin and isospin, although / = J = 1/2 is rather plausible. 
\ Its parity is also not known. Its width (< 10 MeV) is exceptionally narrow as for a hadron resonance located at 
^ . 110 MeV above the nK^ threshold. There have been many theoretical analyses on the pentaquark before and after 
the discovery of 0+ on the basis of the skyrme model, quark models, QCD sum rules, lattice QCD simulations and 
so on, which are summarized in recent reviews 0, . 
QQ The question we address in this paper is the origin of the narrow width of the pentaquark which is the most peculiar 

feature of this new resonance. Since there is no known selection rule from symmetry to make the width naturally 
small, the narrow width should have some dynamical origin. In the past, there have been several attempts to explain 
\ the narrow width from combinatorial suppression from the spin-flavor and color-factors or from the special spatial 
structure due to diquark correlations. However, no definite conclusion has been reached yet. 
rS I In this paper, we propose a possible dynamical mechanism for the long life time of the pentaquark by treating 
JL its hadronic decay on the basis of the constituent quark model and molecular dynamics. Suppose that the colors 
of the constituent quarks are "well-mixed" inside the pentaquark (the precise meaning of the mixing will be given 
(~| , later), then it takes a long time for the quarks to rearrange their colors, flavors, spins, and spatial positions into two 
color-white clusters, i.e. the nucleon N and the kaon K . The dynamical time scale of this rearrangement is governed 
by the strong interactions among quarks and is not simply related to the "distance" in color-flavor-spin-space between 
the pentaquark and the NK state. Therefore, it is not trivial from the outset that the life time is long. 
5-H ' The color molecular dynamics (CMD) simulations originally developed in ref. 6] gives us a suitable framework to 
study the dynamical decay process of the pentaquark. The CMD is a quantum molecular dynamics for constituent 
quarks, in which single quark wave function is parameterized by a Gaussian wave pocket in coordinate space and by a 
color coherent state in the SU (3) color space. Time-dependent dynamics of the multi-quark system is then governed 
by the Hamiltonian commonly used in the standard constituent quark models 0| • The clusterings and decays of the 
multi-quarks are easily treated by this approach. As a first attempt, we will treat the most essential part of the decay 
process (color and spatial coordinates) and neglect spin, flavors and antisymmetrization. 
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II. BASIC FORMULATIONS OF CMD 

We express a total wave function of the system ^ as a direct product of single-particle quark wave- functions: 
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Here i specifies the quarks and anti-quarks. N is the total number of quarks-|-anti-quarks in the system {N = 2, 3 
and 5 for the kaon, the nucleon, and the pentaquark, respectively). <j)i is a Gaussian wave packet centered around Ri 
with momentum and a fixed width L, and Xi is a coherent state in the color SU(3) space parameterized by four 
angles, ai, l3i,9i,ipi . 

Time evolution of the system is given by solving the equations of motion for {Ri, Pi,ai, (3i, 9i, Lpi\ obtained from 
the time-dependent variational principle {5L = 0) on the expectation value for the Lagrangian: 



L = imh H\-^) 
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where H — {'^\H\'^). The Hamiltonian of the system is given as: 
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where t"" = A°/2 for quarks and i° = — A*"/2 for anti-quarks with A° being Gell-Mann matrices. Vc consists of 
one-gluon exchange and confinement terms with an infrared cutoff at r = S.Ofm Typical values of the parameters 
in the quark model for baryons and mesons read niu^d — 300 MeV, = 500 MeV (the constituent-quark mass), 
as = 1.25 (the QCD fine structure constant), and K — G.75GeV/fm (the string tension). We take (size of the 
quark wave packet) to be 0.49 fm, so that the root-mean-square (rms) radius of the quark becomes 0.6 fm which also 
corresponds to the rms radius of the quark core in the ground state of baryons. 

Color-dependent potential Vc{r) alone does not lead to the mass difference between 6+ and nK^ state. Since 
the 110 MeV mass difference is essential for the pentaquark to decay, we introduce V}\f (r) to reproduce this mass 
difference. We take color independent potential with a combination of the attractive scalar-type and repulsive vector- 
type Yukawa potentials P|. We assume that Vj\/f(r) does not act on the s quark. lOj The effective size of the wave 
packet in the matrix element of Vm is chosen to be 0.6 fm to reproduce the mass splitting between 6+ and NK. 



A. Criterion for color- whiteness 

To make a color-white hadronfll| from arbitrary initial configurations, we use a cooling technique in color space 
. Whether the system of five quarks (four quarks labeled by i = 1 ~ 4 and one anti-quark labeled hy i — b) become 
color-white or not is decided by solving the following criterion; 

4 

E(^'I^"I^»)-(X5|A™|X5)=0 (a = l,...,8). (7) 

1=1 

Criterion of color whiteness for the qq and qqq systems are similarly defined. The color-white pentaquark satisfying 
eq.Q is not necessarily composed of two color- white subclusters. This is why dynamical rearrangement of internal 
colors is required for the pentaquark to decay into the NK state. 
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III. DECAY OF THE PENTAQUARK 

A. Color mixing rate (a) and effective spatial distance (D) 

The decay of the pentaquark to the NK state is a dynamical process in which all the coordinate values (10 for 
each quark and 50 for the five quarks) defined in the previous section change in time. Although CMD simulations 
trace time-development of all these coordinates, it is convenient to define some particular combinations which can 
qualitatively characterize the decay process, a and D defined below are such key parameters. 

a is a measure how well the colors are mixed among five quarks. We take notice of an anti-quark in the system and 
calculate a distance between the anti-quark and one of the other quarks in color-space. The minimum value within 
four possible distances is called a. The definition of a in color space can be seen from the following formula: 

« = ^^lilfl - (X5|A*1X5)]'| . (8) 

For a = 0, not only the five quarks are color-white as a whole, but also there are color-white qq and qqq subclusters. As 
a increases, the mixture of color in the five quarks increases. The system cannot be separated into two white clusters 
any more for large a. Note that this parameter is defined only in the color space. Therefore, even if subclusters 
formed at small a are color-while, they do not have to be clusters in coordinate space. 

Next, we define an effective distance D which is a measure the rate of clustering of the five quarks in the coordinate 
space. We consider the anti-quark. A quark closest to the anti-quark is supposed to form a subcluster with the 
anti-quark in the coordinate space. Then, we define a center of mass coordinate as Rqq. Together with the center of 
mass coordinate of the remaining three quarks Rqqq, D is defined by 

D=\Rqq-B.qqq\. (9) 

For small D, the five quarks are in one unit, while for large Z?, they split into two spatially separated subclusters. 
Since D is defined only in the coordinate space, each subcluster is not necessarily color-white. 

By using these parameters, the pentaquark 0+ is characterized as a state with large a and small D, while the NK 
scattering state is characterized by small a and large D. Both parameters are time-dependent and their initial values 
at i = are defined as 

ainit = a(t = 0), Anit = 0). (10) 

We simulate the time development of the five quarks from initial conditions with various color and spatial configu- 
rations (ofinit and -Dinit)- Then we estimate the lifetime of the system until its decay into the NK state. 

B. Correlation between the color mixing and the life time 

Let us first investigate the correlation between the initial color mixing rate and the lifetime. We start with initial 
variables {R^, a^, 0^, (pi} randomly chosen under the constraint that the five quarks are color white as a whole. 
We choose the initial momentum = 0. Then, the equation of motion for 50 (5 x 10) coordinates are solved until 
the system decays into color-white qq {K) + color-white qqq (TV) system. The life time is defined when a becomes 
sufficiently small and D becomes sufficiently large so that the system becomes the NK state. Actual conditions we 
use are 

a < 0.05 and D > 3.0fm. (11) 

Figure^a) shows a relation between ainit and the lifetime T. The solid point for a given T is obtained by averaging 
over randomly distributed values of ainit- The error-bars show the la variance of ainit- This figure clearly shows a 
positive correlation between T and ainit- The result can be easily understood: it takes more time for the state with 
large color mixing to rearrange their colors to white and white, which decays into the NK final state. 

C. Correlation between spatial clustering and the life time 

Next we investigate the relation between the initial spatial distance Dinit and the lifetime T. Figure ^i^b) shows 
the relation between Z?init and the lifetime T. The samples are the same as in the case of Sec. 3. 2. The solid point for 
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FIG. 1: (a) Relation between the initial color mixing rate ainit and the lifetime T. (b) Relation between the initial spatial 
distance Dinit and the lifetime T. 



a given T is obtained by averaging over randomly distributed values of I?init- The error-bars show the la variance 
of the I?init- The figure shows a negative correlation between T and Mnit- The result can be easily understood: it 
takes more time for the state with small spatial distance between subclusters to rearrange their coordinates to the 
NK final state. 

If we choose ainit — 0.5 and I?init = as a typical parameter set for the pentaquark state 0+, we have the life 
time (decay width) of about 100 fm/c (2 MeV). It is too hasty at the moment to compare this number with the 
experimental upper bound of 10 McV, but the result is suggestive. 



D. Potential surface as a function of a and D 



To study the effect of color mixing rate a and the effective spatial distance D on the dynamical decay process in 
more details, we calculate an effective potential V between the qq and qqq subclusters as functions of a and D. 

Shown in Fig. EJa) is the result of such calculation using the potential part of the Hamiltonian in eq.lO- In practice, 
V{a,D) is calculated as follows: first, we prepare a qq state with quark and anti-quark sitting on top of each other 
in coordinate space, and qqq is prepared in the same way. Then, we calculate a, D and the expectation value of the 
potential energy of this state. Then, this process is repeated to cover the variety of points in two-dimensional a — D 
plane. With this procedure, V is adjusted to be zero for the NK final state at large D and small a. Also V^=110 
MeV for small D and large a corresponds to the ideal excitation energy of the pentaquark state 0+. 

In general, V increases as D increases due to the effect of color confinement potential. Namely, it takes more and 
more energy when one tries to separate the five quarks into color non-white subclusters in the coordinate space. The 
only exception is a = where the potential energy decreases as D increases. This is because there is no resistance 
from the confinement force in this case. Note also that the potential surface is fiat in the a direction for small D. 
(One can prove that V is exactly a- independent for D = Q). This implies that the color-dependent potential become 
important only when quarks are separated in space. 

Now, suppose that the pentaquark state 0+ is located at small D and large a region as indicated by the solid circle 
in Fig. [2Ia). For this state to decay into the NK state indicated by the open circle in Fig. |2Ia), it has to rearrange 
color to find a narrow channel near a = 0. This takes a long time since the potential surface is flat in a direction and 
thus the system goes back and forth before reaching the channel. Once it reaches the region around a = 0, it quickly 
decays into the NK state along the narrow channel. Namely, the flat potential along the a-axis near D = and the 
narrow channel along the Z)-axis near a = are two essential sources of the long life time of the pentaquark. 

Figure|3b) shows the actual path in the simulation of the pentaquark state decaying into the NK state. The initial 
conditions are taken to be ainit = 0.4 and Dinit = 0.25 fm, and it takes 25 fm/c to decay. Wandering in color space of 
the pentaquark before the decay can be seen explicitly from this figure. 
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FIG. 2: (a) An effective potential V on the a — D plain, (b) Actual motion of the pentaquark decaying into the NK state is 
superimposed on the top view of the effective potential surface. 



IV. SUMMARY 

In order to study the narrow decay width of the pentaquark suggested by experiments, we have carried out color 
molecular dynamics (CMD) simulation for five quarks. In CMD, the spatial and color parts of the wave function for 
each quark is expressed as the Gaussian wave packet and SU{3) coherent state respectively. An advantage of this 
approach is that one can trace the rearrangement process of spatial and color coordinates as a function of time during 
the decay. 

To characterize the essential part of the decay process, we have introduced two key parameters, the color mixing 
ratio a and the effective spatial distance D. The results of the simulation show that there is a positive (negative) 
correlation between a{t = 0) {D{t = 0)) and the life time T of the five-quark state. T can reach to even 100 fm/c if 
a{t = 0) is enough large and D{t = 0) is enough small. Narrow channel in the effective potential surface V{a, D) is 
found to be the physical origin to cause the long life time. The pentaquark wanders around the potential surface. 

In the present paper, we have not considered the spin and fiavor of the quarks, and hence the spin and flavor 
dependent interactions in the Hamiltonian. The antisymmetrization of the quark wave functions are also neglected 
for simplicity. Detailed study of the decay process taking into account these ingredients is an urgent problem to be 
examined. 

In this paper, we have simply assumed that the pentaquark 0+ corresponds to a state with a large a and small 
D. Whether one can indeed reach such parameter region by the process N + K ^ 8+ N + K can be studied by 
the scattering simulations in CMD, which will be reported elsewhere @. Applications of the present study to the 
decay processes of other possible narrow resonances such as charmed pentaquark are also the interesting problems to 
be investigated. 
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